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ABSTRACT 

If gamma-ray burst (GRB) emission is strongly collimated then GRBs occur throughout the 
Universe at a rate much higher than is detected. Since the emission from the optical afterglow 
is thought to be more isotropic than the gamma-ray emission, it has been hypothesized that a 
search for orphan afterglows (those without the triggering GRB) would allow strong constraints 
to be placed on the degree of GRB collimation. We show here that, within the context of leading 
models of GRB jet evolution, measurement of the GRB beaming angle 9 Jct using optical orphan 
searches is extremely difficult, perhaps impossible in practice. This is because in the leading 
model of GRB jets, the effective afterglow beaming angle scales with the jet angle, f2 op t oc Oj ot 
for small angles, and so the ratio of detected orphan afterglows to GRBs is independent of the jet 
opening angle. Thus, the number of expected afterglow detections is the same for moderate jet 
| angles (~ 20°) as for arbitrarily small jet angles (<C 0.1°). For nearly isotropic GRB geometry, or 

for radio afterglow searches in which the jet has become non-relativistic, the ratio of afterglows 
to GRBs may give information on collimation. However, using a simple model we estimate the 
expected number of orphan detections in current supernova surveys, and find this number to be 
less than one, for all jet opening angles. Even for future supernova surveys, the small detection 
. rate and lack of dependence on collimation angle appear to ruin the prospects of determining 

Q_i' GRB beaming by this method. Radio searches may provide the best hope to find the missing 

q ■ orphans. 

CO 

1. Introduction 

Gamma-ray bursts (GRBs) are observed at a rate of about one per day and are in some cases accompa- 
nied by optical afterglows that are useful in determining properties of the burst such as its redshift. Current 
models that attempt to explain the enormous apparent energy release, the multiband spectra and the tem- 
poral behavior of afterglow lightcurves typically invoke a highly collimated relativistic jet which is beamed 
towards the Earth. The lightcurve from a highly collimated jet is expected initially to decay as a low power 
of t (e.g. flux ~ i -1 - 1 ) and then break to a steeper slope (e.g. flux ~ t~ 2A ) at late times (Piran 1999). A 
number of GRB afterglows show evidence for just this behavior (Sari, Piran & Halpern 1999). This has led 
Frail et al. (2001) to conclude that the GRB emission is collimated into opening angles between 3° and 25° 
with ~ 4° as an average. 

One important feature of the jet model of GRBs is that the vast majority of GRBs are beamed in 
directions away from the Earth and therefore not observed in gamma-rays. However, the afterglow emission, 
occurring at much later time and at longer wavelength should be beamed into a larger fraction of the sky due 
to the decay of the jet Lorentz factor with time. There should therefore be afterglows that are observable 
when the associated GRB is not (Rhoads 1997; Perna & Loeb 1998; Meszaros, Rees & Wijers 1999; Rhoads 
2000). As first suggested by Rhoads (1997), the ratio of the "orphan" afterglow rate to the GRB rate might 
allow measurement of the GRB collimation angle 0j e t, via the equation O opt / S^jct = N gv ^/ N opt , where f2j e t 
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is the solid angle subtended by the jet, fi op t is the angle into which the afterglow is beamed, and N opt is the 
efficiency corrected event rate for optical afterglows. Thus, by counting the rate at which orphan afterglows 
are detected, we could measure the average GRB jet angle, and independently check other methods (Frail 
et al. 2001) of measuring GRB collimation. 

Rhoads' method has already been applied several times to extant searches for type la supernovae and 
other optical transients. Rhoads argued that the lack of afterglow detections by the Supernova Cosmology 
Project and the High-z Supernova Search roughly limited Q opt /tt- )ct < 100, based upon a lack of detections 
over 2 years' exposure. The continued non-detection of orphans since that time has been used by several 
groups to set stronger collimation limits. For example, Rees (1999) found O opt /fij e t < 20, and Rhoads (2001) 
quotes a transient search by Shaefer et al. to claim f2 pt/fijct *C 100. If optical afterglow emission is only 
mildly beamed (e.g. f2 op t « 1), then this is marginally inconsistent with the collimation factors (~ 500) 
claimed by Frail et al. (2001), and could render untenable most solar mass progenitor models, given the 
inefficiency of converting rest mass into gamma-ray emission (Kumar 1999). 

In this paper, however, we argue that current supernova searches place no constraints upon GRB 
collimation. One reason for this is that present surveys are too shallow, narrow, and infrequent to detect 
significant off-axis emission. Using a simplified model of GRB afterglow emission and survey detection 
efficiency, we show that present supernova surveys are generically expected to detect (far) less than one 
orphan afterglow. A deeper reason for the impotence of present surveys is the fact that, within the context 
of the leading GRB afterglow model (Sari, Piran & Halpern 1999; Rhoads 1999), the effective afterglow 
beaming angle scales with the GRB jet angle, opt oc 0j ct . Because of this scaling, the ratio O opt /Oj et is 
unchanged even as the GRB jet angle is varied from moderate angles (0j ct ~ 20°) to arbitrarily small angles 
(#jet -C 0.1°), and so the ratio of detected orphan afterglows to GRBs cannot be used to measure 6>j ot . 

We derive this scaling below using a simplified model of GRB jet emission, however we can understand 
this result in more general terms. From relativistic kinematics, we know that the jet Lorentz factor 7, 
beaming angle 0, and observer time t are related by simple power laws, and as mentioned above the observed 
flux / is a simple power law in time. Thus, the beaming angle is a (scale free) power law in 7, t and /. In 
the limit 7> 1 the only angular scale in the problem becomes 0j et , so we expect the effective beaming angle 
to scale as the jet angle, opt oc 6*j e t- This should break down in the nonrelativistic limit 7 ~ 1. 

This scaling of op t with 0j et severely hampers the ability of optical orphan searches to measure or 
constrain GRB collimation. In fact, we argue in this paper that measurement of the GRB jet angle may 
be impossible using either present or future optical supernova surveys. The plan of the paper is as follows. 
In § 2 we describe our simple model of emission from a GRB afterglow. In § 3 we estimate the afterglow 
detection rate using this simple model. In § 4 we discuss the results of our calculations, and in § 5 we give 
our conclusions, along with some discussion and caveats. 



2. A simple model of optical afterglow emission 

GRB-alerted afterglows ("on-axis afterglows") arc now well studied. In most cases, their time depen- 
dence can be characterized as a broken power-law, / ~ i~ Q , with a power-law break from a « 1.1 toasi 2.4 
occurring at late times when the relativistic beaming angle 7 -1 exceeds the jet opening angle 6*j c t- Here 7 is 
the bulk Lorentz factor of the relativistic jet. This break occurs at a time (Sari, Piran & Halpern 1999) 

i brcak » 6.2(£ 52 /ni) 1/3 (0jct/O.l) 8/3 hr, (1) 
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where E 52 is the inferred energy of the ejecta assuming isotropic expansion, in units of 10 52 ergs, and n\ is 
the density of the surrounding ISM in cm~ 3 . GRB optical emission has a wide range in apparent brightness, 
ranging from 9th magnitude for the prompt optical emission associated with GRB 990123, down to the 
limiting magnitude R ~ 23 — 24 of follow-up searches for afterglows. Fynbo et al. (2001) have argued that 
most (~ 70%) afterglows are even fainter than this, and have therefore eluded detection. 

Observers of orphan afterglows, on the other hand, would see off-axis lightcurves, which differ from 
on-axis lightcurves. The behavior of off-axis emission is governed by relativistic effects, as well as the lateral 
spreading and internal structure of the jet. We will first consider the effects of relativistic kinematics, and 
then consider lateral spreading and jet internal structure. Since the internal structure and time evolution 
of GRB jets is quite uncertain we will consider throughout two limiting cases which we hope will bracket 
reality. 

As a simple model, consider a radiating plasma moving relativistically along the jet axis, with bulk 
Lorentz factor 7 relative to us, and radiating isotropically in its own frame. The measured power, as a 
function of angle 9 with respect to the jet axis in our frame, is (Rybicki & Lightman 1979) 

^ = I (2) 

dQ 7 4 (1 -/?cos6>) 4 V ; 

where F is the (isotropic) power per unit solid angle in the plasma frame, and = (1 — 7~ 2 ) 1//2 as usual. 
Let us assume that the Lorentz factor 7 and rest frame flux F evolve with time to give an observed on-axis 
(9 = 0) lightcurve f(t) matching the lightcurves observed for optical GRB afterglows, i.e. f(t) ~ t~ a . For 
early times before the break a = l-1.2«l.l while at late times after the break a w 2.4. Since 7 oc with 
ji = 3/8 before the break and fi — 1/2 after the break (Sari, Piran & Halpern 1999), this means / oc 7"/'*. 
This then requires F oc "f a ^(dP/dfl(9 = 0)) _1 . With our assumption of rest frame isotropy, we obtain the 
flux seen by off- axis observers as a function of angle 9 and Lorentz factor 7, 

/(t) -Mi=r(T^=»)" (3) 

Here /break is the flux seen by an on-axis observer at the break time ibroak, and we have generalized the 
exponent 4 to the parameter m for later use. This function with m = 4 is plotted as a function oft for several 
values of 9 in Figure 1. Several features follow from equation 3 and Figure 1. First, the lightcurve seen at 
angle 9 peaks at a time i pca k which can be estimated by setting 7 « 9~ x . Next, the flux rises to the peak 
roughly as t 2m ^ a , and decays after peak as t~ a . Also, the flux is nearly independent of angle for 9 <C 7 -1 ; 
that is the various lightcurves all match together at late times when the viewing angle becomes small relative 
to the beaming angle. Finally, when 7 < 6* _1 , an off-axis viewer at 9 > 9 ]ct should see comparable flux to an 
on-axis observer inside 6*j et , and the off-axis afterglow lightcurve should decay with roughly the same slope 
as the on-axis afterglow. 

Note that off- axis (9 > 9j ct ) afterglows do not peak until after the on-axis power-law break, which 
typically occurs hours to days after the GRB (Frail et al. 2001). This implies that orphan afterglows, even 
at peak brightness, should be far fainter than their GRB-triggered siblings, and thus should be much more 
difficult to detect. We can roughly estimate how much fainter as follows. After the break, the afterglow flux 
decays as t~ 2A , while the Lorentz factor decays as 7 oc i~ 1//2 (Sari, Piran & Halpern 1999), so the flux scales 
as 7 4 8 . For a reasonable estimate of the jet angle, we would like to detect orphans at large viewing angles, 
say 9 w 0.2 — 0.3. This means detecting emission at times when the Lorentz factor has decayed to 7 « 3 — 4. 
Taking 7 « (4 ) -1 ~ 15 at the break, we see that the flux is 7 or 8 magnitudes fainter by the time 7 « 3 — 4 
than at the break. 
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This immediately shows why current supernova searches are not expected to detect orphan afterglows. 
Typical limiting magnitudes for such searches are near R w 23, which Fynbo et al. (2001) claim is too 
shallow to detect most on-axis afterglows, let alone far fainter off-axis orphan afterglows. This implies that 
detectable orphan afterglows are mostly contained in a cone of opening angle not much larger than 0j et , so we 
do not expect many more orphan afterglows to appear in SN surveys than on-axis afterglows. Additionally, 
we show in the next section that the relatively infrequent monitoring characteristic of current SN surveys 
(e.g. i mon = 21 days) causes most bright orphans to be missed. 

We shall quantify these broad statements in the next section with a calculation of the detection rates 
expected in surveys, using this simple model of off-axis afterglow emission. Before proceeding, we pause to 
note that we have, so far, neglected the effects of finite jet width and jet spreading. Both of these should 
lead to greater off- axis emission than indicated by equation 3. To obtain an upper limit to the effects of the 
uncertain jet evolution (and thus bracket the range of possible fluxes seen by off-axis observers), we consider 
the following simple picture. The effect of finite jet width on the off-axis lightcurve is to point some of the 
emitting material more closely to the observer's line of sight, and to point other emitting material further 
from the observer's line of sight. A proper calculation of this effect would involve specifying the emission 
profile in the jet and integrating over the jet face (Woods & Loeb 1999), but in the spirit of setting upper 
bounds on the detection rate, we can get a sense of the effect by simpler means. Consider the emission from 
the part of the jet face most closely pointing to a fiducial off-axis observer at 9. This material is moving 
at an angle 9 — 6*j ct from the observer's line of sight, so we can approximate the flux from this part of the 
jet face by replacing 9 with 9 — 9- iCt in equation 3. We immediately see that the asymptotic scalings remain 
unchanged (i.e. t 4 ~ a and t~ a ), but that the peak time changes, to w 9 — 9 }ct . This makes sense, because 
we roughly expect the flux to beam into an angle ~ 9j Ct + 7 _1 - This of course overestimates the total flux 
and the true visible angle, because emission from the rest of the jet face is less blueshifted and therefore less 
bright. Also note that the emission from other regions of the jet face should peak at later times, so that we 
expect a broader lightcurve peak for off axis viewers than on axis viewers. This broadening of the peak is 
seen in numerical simulations of relativistic jet hydrodynamics and emission (Granot et al. 2001). 

An upper limit on the effects of lateral spreading of the jet on the off-axis lightcurve may be obtained by 
noting that if in the jet comoving rest frame the emitting material is expanding laterally at speed v, then the 
apparent angular size of the jet face becomes 6j et = #j e t,o + (^/c)7 _1 - Then our previous discussion applies, 
merely replacing the constant jet angle with this spreading jet size. To bracket possible jet morphologies and 
evolution, for a lower bound on off-axis emission we use equation 3, which has all emission coming from the 
jet center and no lateral spreading, and for an upper bound on off-axis emission we use equation 3 with 9 
replaced by 9 — 9- iCt — 7 _1 , which has all emission coming from the nearest jet edge and lateral jet spreading 
at the speed of light. 

3. Event detection rate 

Observing programs that aim to measure GRB collimation using orphans will measure the rate of orphan 
afterglow detection and compare this to the rate of GRB-triggered (on-axis) afterglow occurrence. Optical 
searches for orphan afterglows naturally piggyback on top of Type la supernova searches, which repeatedly 
monitor fixed regions of the sky and use differential photometry to detect transient optical sources. Ideally 
one would like to detect every afterglow occurring in the observed region, but some afterglows will be missed 
either because they are too faint or because the monitoring is not frequent enough to catch them. Note, 
however, that if one can accurately estimate the efficiency at which afterglows are found, then orphan searches 



- 5 - 



may be useful even if many afterglows are missed. 

In this section, we estimate the number of afterglows that an observing program would detect and the 
properties (such as viewing angle from the center of the jet) of those detected afterglows. We can estimate 
these quantities using the simple model of off-axis emission described in the previous section. First, we 
give a simple analytic calculation of the detection rate, appropriate in the limit of frequent or continuous 
monitoring, which captures the effects of jet behavior. To determine the effects of infrequent monitoring, we 
then perform a Monte Carlo simulation of a simplified observing program. 



3.1. Analytic estimate 



Consider a perfect orphan afterglow search capable of detecting any orphan afterglow above limiting 
flux /i im . We would like to determine which viewing angles 9 see apparent fluxes above threshold, / > /ii m , 
as a function of time t. Using equation 3 and 7 = #j e t _ (*Abreak)~' 1 > and taking 7 >few, we see that / > /u m 
for viewing angles 6 < # max , where 6* max is given by 
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Recall that before the break, /j, = 3/8 and a « 1.1, while after ibrcak we have [i = 1/2 and a w 2.4. To 
compute the effective afterglow beaming angle, we maximize equation 4 with respect to time. Let us call 
t max the time at which the visible angle is maximized. Depending upon whether or not i max coincides with 
the break time, the # max expression becomes 
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This discussion has neglected both the finite width of the jet and the possibility of lateral spreading of 
the jet. Including these effects as discussed in § 2 gives 



0max(z) — #jct + ^1 + — ^ #max,0 



(6) 



where 6> maXi0 is defined as in equation 5. We emphasize that this is an upper limit to the visible angle, 
and that a more careful treatment would give a # max somewhere in between that given by equation 5 and 
equation 6. Our previous expression for # maXj o should underestimate off- axis emission, and so these two cases 
should bracket all possible jet emission. 

In order to proceed, we next need to define the GRB properties. There are large uncertainties here, 
and we expect these to introduce substantial uncertainties in our estimates of the event rates. However, we 
argue that this will not affect our main conclusions. 

First we follow Wijers et al. (1998) and Woods & Loeb (1998) in hypothesizing that the comoving rate 
at which GRBs occur is proportional to the star formation rate, -Rsfr- For simplicity we approximate the 
star formation rate, as measured by Steidel et al. (1999), as 



logio r sfr{z) 



Cz 
C 



z < 1 

1< z < 10, 



(7) 
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with an arbitrary cutoff for z > 10. For the VLm = 1, Hq = 50 km/s/Mpc Einstein-de Sitter universe that 
Steidel et al. use, C « 1 appears adequate, and correcting this to the fiat = 0.3, Ha — 0.7, H a = 70 
km/s/Mpc cosmology that we adopt gives instead C « 0.75. We absorb the overall normalization of i?sFR 
into our total event rate, which we normalize so that the overall rate of GRBs is <~ 666/yr over the full sky 
(Paciesas et al. 1999). 1 

Next we must make an assumption concerning the luminosity function of GRB afterglows. For simplicity 
we will assume that afterglows are standard candles with absolute magnitude M = — 25 at the time of the 
lightcurve break. This number is chosen so that ~ 70% of well-followed GRBs do not show afterglows 
because the afterglows are too faint at the time of observation, which we assume to be near the break time. 
Of course, it is known that the dozen or so afterglows with measured redshifts are widely distributed in 
absolute magnitude, but M = —25 is a typical value. Since the effective beaming angle depends strongly on 
the apparent luminosity function (c.f. equation 5), uncertainty here will lead to substantial uncertainty in 
the event rate, as we discuss below. 

One additional concern is the k-correction. At high frequencies, the afterglow spectrum is expected to 
behave as F v ~ v~ a / 2 (Sari, Piran & Halpern 1999), so as a function of redshift z, /broak( 2 ) = /break(O) x 

[1 + Z]- Q / 2 . 

Finally, we assume that all GRBs have the same jet opening angle 9j ct , the same power law indices a 
and /« specified in § 2, and the same rest- frame break time tbrcak ~ 1 day. Note that here by rest-frame, we 
mean as viewed by an observer at the same redshift, so that there is a (1 + z) correction. We do not mean 
the plasma rest frame time, for which there is an additional 7 2 correction relative to observer time. Clearly 
there is a wide spread in the observed break's and power-law indices a for previous GRBs (Frail et al. 2001), 
and this will lead to uncertainties of order a few in our final numerical results. 

We can now compute the expected rate of orphan afterglow detections. Let us write the effective 
beaming (solid) angle as f2&(z) = Att(\ — cos(9 max (z)). With our assumption of standard candle afterglows, 
we need only specify the redshift distribution to obtain the afterglow luminosity function. Assuming as 
above that the gamma-ray burst rate is proportional to the star formation rate, we may write the observed 
afterglow redshift distribution as 

dz dz 1 + z 

where dV co /dz is the comoving volume element, z is the GRB redshift, and A is a normalization constant. 
The k-correction is incorporated as above and we set A as above by assuming that all on axis gamma-ray 
bursts are bright enough to be detected; that is we set max = 8j Ct in equation 8 and integrate over redshift 
giving N GKB = A ■ AO jet J °° dz{dV co / dz)R SFK / + z) = 666 yr _1 over the full sky. 

The calculation described above elucidates the overall behavior, but neglects effects which may affect 
numerical predictions for the overall detection rate. For example, we have extrapolated the scaling formulas 
for 7 to the nonrelativistic regime where the scalings break down; fortunately we will show that the regions 
of low 7 make a small contribution to the integral for all but the deepest surveys. We have also assumed 
continuous monitoring, which is certainly not the case in current SNe searches. In order to take into account 
the effect of imperfect monitoring and to check the above analytic calculation we ran a Monte Carlo simulation 



1 Note the number 666/yr includes both long and short GRBs, while afterglows have been detected only for long bursts. We 
conservatively assume here, however, that all bursts give rise to afterglows; this will soon be testable with experiments such as 
HETE-II or Swift. 
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of a simplified observing program. 

3.2. Monte Carlo simulation 

We define an observing program by its limiting magnitude mu m (and corresponding limiting flux /n m ), 
by the time between exposures of a given field t mon , and by the total survey exposure (number of deg 2 yr). 
In an actual survey additional parameters and various experimental efficiencies would also be taken into 
account but these are not crucial for our considerations here. The observing program simulation is run for a 
given value of #j e t by distributing GRBs randomly throughout the Universe, proportionally to the product 
of comoving volume and star formation rate, and at a random viewing angle to the observer. The lightcurve 
is started at a random time and observed every i mon days. The afterglow is considered to be detected if any 
of the observations find a flux above the limiting flux of the observing program. The total event rate and 
the distribution of detected events as a function of redshift, observing angle, etc. is thus easily found. 

4. Results 

Now, all the elements are in place to calculate the detection rate and average beaming factor for a given 
survey. Let us first consider a Type la supernova search (Schmidt et al. 1998; Pcrlmutter et al. 1999), with 
limiting magnitude of R s=s 23 and monitoring repeat baseline t mon « 21 days. Adopting an opening angle 
of 0jet = 4° (Frail et al. 2001) our Monte Carlo calculation gives an expected rate of ~ 0.0008 deg~ 2 yr _1 if 
the flux is concentrated at the jet center and there is no lateral spreading, and ~0.014 deg~ 2 yr" 1 if lateral 
spreading is maximal and the flux is concentrated at the jet edge. So, taking an exposure of 10 deg 2 yr, the 
we predict between 0.008 and 0.14 detections from previous Type la supernova surveys. We note that in 
the limit of very small i mon the Monte Carlo results are close to the analytic results, which finds 3-5 times 
more events depending on jet morphology. 

Now, the important point to note is that these numbers, besides being small, are fairly independent of 
8j ct . The reason for this should be clear. The only dependence of equation 8 on 6j ct is the fib(z) factor, which 
was the effective beaming angle as a function of redshift. Since max cx 6>j ct , fi b cx fij et . But since the total 
GRB rate scales like fir* to hold fixed the number of observed GRBs, this means that the expected number 
of detected orphans is also held fixed. So the non-detection of orphan afterglows is not only consistent with 
opening angles of 4°, it is consistent with almost any opening angle. 

Next we consider a more powerful future observing program with limiting magnitude i?n m w 27, and 
frequent enough monitoring that our analytic estimates are accurate. Here the number of predicted afterglow 
detections can rise above unity if many tens of square degrees are monitored. The number of expected events 
is shown in Figure 2 as a function of the opening jet angle #j e t- For opening angles ranging from moderate 
(~ 20°) to arbitrarily small (<C 0.01°) there is essentially no change in the predicted number of detections. 
Note, however, that for large jet angles (> 30°) there can be a significant decrement in the predicted number 
of detected orphans, in part because the relativistic scalings we have employed break down, and in part 
because the emission begins to subtend the full sky. However, it is clear that for the angles being considered 
in the leading jet models, the ratio of GRB detection to afterglow detection will give no information on the 
GRB collimation. 
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5. Discussion and conclusions 

We have shown that within the framework of the leading model for the late time GRB jet evolution, 
current optical surveys cannot constrain GRB beaming. Part of the problem is that current surveys are too 
infrequent, narrow and shallow. 

A more basic problem is that the effective afterglow beaming angle scales as the jet angle for small 
(< 20°) jet angles, and so the ratio of f2 opt /f2j ct cannot be used to measure 9j ct . This conclusion depends on 
the simple scaling laws of flux and Lorentz factor that apply during the relativistic phase of the expansion. 
In order to see emission from large viewing angles, one must wait until the jet has become only moderately 
relativistic and therefore dim. In Figure 3 we plot the distribution of viewing angles for several different 
limiting magnitudes. Only the R = 30 lightcurve has any sensitivity to large viewing angles, and it is still 
dominated by small 9. Finding optical afterglows this faint or fainter will be difficult due to, for example, 
host galaxy domination. 

Absent the capability to detect emission from viewing angles near unity, measuring the GRB jet angle 
using orphans is all but impossible. This is shown in Fig. 2, in which we plot the number of expected 
afterglow detections per year per square degree as a function of jet opening angle, for a limiting magnitude 
of R = 27. As mentioned above, there is essentially no variation in the predicted number of detections as the 
jet solid angle Oj ct is varied by orders of magnitude, until nearly isotropic jets (6*j e t ^ 30°) are reached. Now, 
once the GRB afterglow luminosity function is well measured, as should be possible with future missions 
such as Swift or HETE-II, and once GRB jet dynamics are better understood, then a prediction can be made 
for the expected number of afterglows any given survey can detect. A significant downward departure from 
this expected rate would be hard to reconcile with small jet angles, within the context of accepted models. 
A statistically significant detection of such a deficit would require a very deep survey expected to detect 
hundreds of orphans over a large region of sky. 

To establish our conclusion of the difficulty of using optical afterglow searches to constrain small beaming 
angles, we have presented two sets of detection rates in an attempt to bracket the effects of uncertain jet 
morphology and evolution. We caution that our upper bound is probably not realistic because it assumes 
that the jet emits all of its flux as near as possible to the direction of the observer. In addition, for both the 
upper and lower bounds we have neglected the frequency redshift of photons received by off axis observers 
relative to on-axis observers. This should be taken into account for observers interested in the narrow 
frequency range within the optical band. For a narrow pencil beam moving along the 8 — axis the relative 
redshift of photons received by an observer at 8 b s is v{9 o b s )/v{0 = 0) = (1 — /3cos# fc s )/(l — (3). If the 
afterglow spectrum behaves as F v oc v~ a / 2 (Sari, Piran & Halpern 1999), then the redshift effect drives the 
exponent in equation 3, which we previously found to be m = 4, to m = (4 + a/2). Replacing m = 4 in 
our equations with m = 4 + a/2 slightly diminishes the expected number of detections by < 30%. Were an 
orphan afterglow search to be performed, these corrections, as well as a detailed study of the jet evolution 
and emission, would be required to calculate a more accurate estimate of the expected rate. 

In our calculations we have assumed particular values for the parameters describing GRB afterglows. 
These parameters are the power law exponents a and ji, the absolute magnitude of the on-axis afterglow at 
break, and the break time ibreak- Uncertainties in these parameters will introduce uncertainties in the overall 
rates. For example, changing a from 2.4 to 1.8 increases the detection rate by a factor of ~ 3. Similarly, 
changing the absolute magnitude at break from M = — 25 to M — —29 (an increase in luminosity by a 
factor of 40), increases the overall detection rate by a factor of 4-5, as expected from equation 5. Changing 
ibreak does not affect our calculations for continuous monitoring, but can change the detection efficiency 
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depending upon the monitoring time. In the limit of small efficiency the rate should scale as ibrcak- Most 
importantly, the constancy of the calculated detection rate with respect to the jet opening angle remains 
valid for jet < 20°. 

One caveat to our estimate for the orphan afterglow detection rate and to our constraint on the observed 
ratio of orphans to GRBs is the possibility of processes that lead to weakly beamed optical emission while 
remaining consistent with the good agreement between the expected and observed power law decay and 
break of the late time afterglow. For example the GRB central engine might emit, coincident with an ultra- 
relativistic jet, a moderately relativistic spherical shell that leads to isotropic optical emission. This might 
occur in central engine models where rotation leads to small baryon contamination and ultra-relativistic 
flow along the rotation axis, but also to "dirtier" , moderately relativistic flow away from the rotation axis 
(e.g. in the collapsar model, MacFadyen & Woosley 1999). Such an isotropic component, if bright enough 
to be detected off-axis but faint enough to avoid on-axis detection, would clearly confound our analysis 
and allow direct measurement of GRB collimation. On the other hand, it seems unlikely that we have 
underestimated the off axis optical emission from the jet itself, because this estimate derives principally from 
simple relativistic kinematics. Similarly, the upper limit to the degree of angular spreading in the jet comes 
from the maximal assumption that the lateral expansion velocity is c in the jet comoving rest frame. Recent 
numerical hydrodynamics simulations (Granot et al. 2001) of the relativistic jet dynamics appear to indicate 
that the afterglow emission is indeed more tightly beamed in the forward direction than we have assumed, 
suggesting that our calculations neglecting spreading are more realistic than our calculations including lateral 
spreading. 

Finally, we note that radio afterglows are not subject to the analysis presented here. They do not 
decay as steeply as their higher frequency counterparts and remain visible until the afterglow emission is 
effectively isotropic. Perna & Loeb (1998) have proposed that a full-sky survey sensitive to 0.1 mJy could 
legitimately test the jet model of GRB emission. It appears that this may be the best way to implement 
Rhoads' suggestion of using burst-less afterglows to measure the GRB jet angle. 
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Fig. 1. — Apparent magnitude of a GRB afterglow 
as a function of time for various viewing angles as 
given by equation 3. From top to bottom the curves 
are for viewing angles of 9 =0.01, 0.05, 0.1, 0.2, and 
0.5 radians, respectively. The GRB is taken to be at 
a redshift of 1, have absolute magnitude M = —25 
at the break time of 1 day in the GRB rest frame 
(corresponding to observed break time of ibrcak = 2 
days). The jet opening angle is 9 ]C t = 0.07 radians, 
and it is assumed that the flux is concentrated at 
the jet center and lateral spreading of the jet is not 
included. Note that the flux peaks near the time 
when 7 ps 6~ l . Before the peak, the flux rises as 
roughly t 2m ^~ a , and after the peak it decays as t~ a . 
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Fig. 2. — Plot of total number of afterglows detected 
per year per square degree as a function of jet an- 
gle 0jet in radians for our simple approximation. A 
limiting magnitude of R — 27 was assumed. The dif- 
ferent curves show the effects of jet width and lateral 
spreading. The top two curves assume emission from 
the near edge of the jet (equation 6) with v = c and 
v = respectively, while the bottom curve assumes 
emission only from the jet center (equation 5). 
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Fig. 3. — Plot of dN/d9, in units of events/yr/deg 2 , 
as a function of viewing angle 9, for an assumed jet 
angle 9 ict = 4° and for our analytic approximation. 
The different curves within each panel show the effect 
of varying the limiting magnitude. From top to bot- 
tom, they are mn m = 30, 27, 24, and 21, respectively. 
The different panels correspond to different assump- 
tions of the effects of jet width and lateral spreading. 
In the first panel, off axis emission is estimated using 
equation 5, that is, it assumed to come only from the 
jet center. In the second panel, off-axis emission is es- 
timated using equation 6 with v — c, corresponding 
to emission from the nearest jet edge and maximal 
spreading. Note that dN/d9 is independent of limit- 
ing magnitude for 9 < 9j ct when we include the finite 
jet width, since on-axis afterglows in our model can 
be arbitrarily bright at early times, and we assume 
continuous monitoring in our analytic approximation. 
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